theta = linspace(0, 2*pi, 1000);
z = exp(1i*theta) - 1;
plot(real(z), imag(z)); 

hold on
z = 2 * (exp(2*1i*theta) - exp(1i*theta))./(3*exp(1i*theta) - 1); 
plot(real(z), imag(z)); 

hold on
z = 12*(exp(3*1i*theta) - exp(2*1i*theta))./(23*exp(2*1i*theta) - 16*exp(1i*theta)+ 5) ;
plot(real(z), imag(z));
hold off

z = 24*(exp(4*1i*theta) - exp(3*1i*theta))./(55*exp(3*1i*theta) - 59*exp(2*1i*theta) + 37*exp(1i*theta) - 9) ; 
plot(real(z), imag(z)); 

axis equal; % 设置坐标轴比例相等
xlabel('Real'); % 设置x轴标签
ylabel('Imag'); % 设置y轴标签
